Distribution of the largest fragment in the Lattice Gas Model 
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The distribution of the largest fragment is studied in different regions of the Lattice Gas model 
phase diagram. We show that first and second order transitions can be clearly distinguished in 
the grancanonical ensemble, while signals typical of a continuous transition are seen inside the 
coexistence region if a mass conservation constraint is applied. Some possible implications of these 
findings for heavy ion multifragmentation experiments are discussed. 
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I. INTRODUCTION 

°' . 

Since the first heavy ion experiments the size of the largest cluster Am detected in multifragmentation events 
ft , has been tentatively associated to an order parameter for the fragmentation phase transition^]; if this is true, we 
should expect for this observable a double humped distribution if the transition is first order 0, while its fluctuations 
Ph ', should obey the first scaling law if the transition is continuous 0. Experimental multifragmentation data show in this 
respect somewhat contradictory evidences. An analysis of 80 A.MeV Au+Au peripheral collisions from the Indra- 
i—{ ' Aladin collaboration |4| reports a bimodal distribution of a variable closely correlated to Am- On the other side the 
functional relationship between the two first moments of Am in central Xe+Cu collisions^! shows a change of slope 
which has been interpreted as a a transition from the A = 1/2 to the A = 1 scalin g la w as expected for a generic 
continuous transition^. From the theoretical point of view, it is well known|f3.M.l8ll^.llOl| that in finite systems many 
different pseudo-critical behaviors can be observed inside the coexistence region of a first order phase transition. In 
particular concerning the order parameter fluctuations, simulations have been performed in the framework of the Ising 
■ Model with Fixed Magnetization (IMFM) in ref.^J- I n this study the distribution of Am was shown to approximately 
obey the first scaling law even at subcritical densities, i.e. in thermodynamic conditions where no continuous transition 
l/^ , takes place. Since the scaling is violated for very large lattices, the observed behavior was interpreted in this paper 
as a finite size effect that prevents to recognize the order of a transition in a small system. An important difference 
subsists though between the theoretical study of ref.^l| and the experimental analysis ref.Q: in the first paper the 
average Am size is varied by increasing the total lattice size, meaning that the existence of a scaling law is tested in 
well defined thermodynamic conditions (a single point in the state variables space). In the experimental case 

it is not possible to freely vary the source size, therefore different regions of (Am) are explored by varying the total 
energy deposited in the fragmenting system. It is not a priori clear how these two very different procedures might be 
related and whether they could be equivalent. 

In this paper we analyze the distribution of Am within the Lattice Gas Model^3- This model is the simplest 
representation of the liquid-gas phase transition; once augmented with the cluster definition through the Coniglio- 
Klein algorithm |13|, it can also be related to a bond and site percolation problem, making this model a paradigm of the 
fragmentation phase transition. This model is isomorphous to the Ising spin model and its thermodynamic properties 
are very precisely known: the Lattice Gas phase diagram contains both first and second order phase transitions and 
basic effects, like conservation laws, which are very relevant to the experimental situation, can be easily implemented. 

In the analysis of the Am distributions we will show that the most important finite size effect is the inequivalence 
between statistical ensembles the observed ambiguities can be coherently interpreted as an effect of conservation 
laws, the distribution of an order parameter being drastically deformed if a constraint is applied on an observable 
which is closely correlated to the order parameter under study. 
Specifically we will show that: 

• in small canonical systems, a first scaling law as a function of the system size can be observed at the critical point 
but also for subcritical densities inside the coexistence region. This is in agreement with the findings reported 
in ref. [HI- The difficulty in recognizing the order of the transition is not only due to the finite size effects but 
more important, the order parameter distribution and its scaling properties are deformed by the conservation 
law that in the canonical ensemble acts on the total number of particles A t , which strongly constraints the 
order parameter Am', 
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if the Am size is varied by changing the system temperature at a fixed lattice size, no scaling of the largest 
fragment distribution is observed even if we choose a transformation which passes across the thermodynamic 
critical point; 
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• in this case, the correlation between the average < A m > and the variance a 2 of the largest fragment distribution 
exhibits a rise and fall which is imposed by the conservation law constraint; the double logarithmic derivative 
A' = d log a/d < Am > appears to be a smooth decreasing function of < Am >', even if A' is passing through 
A' = 1 and A' = 1/2 before becoming negative no simple scalings a oc< Am > A can actually be unambiguously 
isolated. 

• however, we show that both the existence of a transition and a conclusion about its order can be infered from 
the quantitative study of the Am fluctuation. 

II. PHASE TRANSITION IN THE LATTICE GAS MODEL 
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FIG. 1: Right part: distributions of the size of the largest cluster in the grancanonical lattice gas model at different temperatures 
T/e, with a 8x8x8 lattice and at [i = 3e. Left side: lattice gas phas diagram from the distributions on the right side. Dashed 
line: locus of the maximal Am fluctuation in the canonical ensemble. 

In our implementation of the lattice gas model [12J the N sites of a cubic lattice are characterized by an occupation 
number m which is defined as rii — 0(1) for a vacancy (particle). Particles occupying nearest neighboring sites interact 
with a constant coupling e. This model can be transformed into an Ising spin problem with a magnetic field through 
the mapping Sj = n, — 1/2. The relative particle density p/po is defined as the number of occupied sites divided by 
the total number of sites and is linked to the magnetization of the Ising model by p/po = m+ 1/2. In addition to this 
interaction part a kinetic energy is introduced. Occupied sites are characterized by a momentum vector. Observables 
expectation values are evaluated in the different ensembles (grancanonical, canonical and microcanonical) sampled 
through standard Metropolis algorithms Jj|. The chemical potential in the grancanonical implementation plays the 
role of the magnetic field h = p—3e in Ising, while the canonical Lattice Gas corresponds to the constant magnetization 
Ising IMFM case with m — p/p — 1/2. 

The phase diagram of the model can be easily evaluated looking at the distribution of the total number of particles 
A t = y\-_i rii in the grancanonical ensemble with a chemical potential p — p c — 3e which corresponds to the Ising 
critical field h = 0. The A t distributions, Pp^(A t ) : are displayed at different temperatures in the right part of figure 
^ The presence of two different ensembles of states (bimodality) is clearly seen for all temperatures T <T c k, 1.22e. 
At the critical chemical potential p c presented in the figure, the probabilities of occurrence of the two solutions are 
exactly identical; if p < p c (p > p c ) the high (low) density peak dominates. For a fixed temperature , the most 
probable A t as a function of p is discontinuous at the transition point p c . 

At the thermodynamic limit, the discontinuity in the most probable A t as a function of p give rise to a discontinuity 
in the associated (A t )(p) equation of state; this implies that the two peaks represent two coexisting phases 0,^^ and 
that the number of particles ( or equivalently the density) is the order parameter of a phase transition which is first 
order up to the critical point T = T c . The phase diagram can be constructed by reporting the forbidden region for 
the most probable density i.e. the locus of the discontinuity in the most probable. This corresponds to the two peaks 
in the bimodal particle number distribution observed at p c . The phase diagram is displayed in the left part of figure 
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^ These findings obtained in a 8x8x8 lattice correspond to the phenomenology of the liquid gas phase transition 
that the model is known to display at the thermodynamic limit. If we increase the lattice size the location of the 
coexistence border will be modified, even if finite size corrections are especially small in this model@. However it is 
clear from figure ^ that (except at the critical point which is a second order point, where the two peaks merge to form 
a single distribution) the first order character of the transition is indisputable even for a linear dimension as small as 
L = 8. 




Temperature/^ 

FIG. 2: Am distributions as a function of temperature for a 8x8x8 lattice in the grancanonical (left), canonical (middle), and 
microcanonical (right) ensemble. In all cases the density is p/po = 1/2. 

Figure|21shows the size of the largest cluster Am as a function of the temperature for the grancanonical, canonical and 
microcanonical ensembles. The obvious correlation between Am and A t implies that for T < T c the Am distribution 
is also double humped in the grancanonical ensemble as explicitly shown in ref.[17|. This means that Am can also be 
taken as an order parameter of the liquid gas phase transition, and looking at its distribution this transition can be 
recognized as first order even for a system constituted of (At) — 256 particles. 



III. CONSERVATION LAWS AND THERMODYNAMICS 



If the constraint of mass conservation is implemented (canonical lattice gas, or equivalently Ising model with fixed 
magnetization) the distributions of Am drastically changefl7j. In the grancanonical ensemble at (p/po) = 1/2, the 
explored microstates essentially populate the coexistence border, while the coexistence region is accessed with a 
negligible probability (see figure these highly improbable of the grancanonical distributions are conversely the only 
microstates which are allowed by the canonical constraint at the value p/po = 1/2 below the transition temperature 
the grand canonical and canonical partitions differ drastically. Because of the mass conservation constraint, the 
bimodality of the A t distribution is obviously lost in the canonical ensemble; as a consequence of the correlation 
between A t and Am, the Am distribution also shows a unique peak (figure^. If we additionally implement a total 
energy conservation constraint (microcanonical ensemble, right part of figure [21 the distributions get still narrower, 
but the qualitative behavior is the same than in the canonical ensemble. The normal behavior of the Am distribution 
at subcritical temperatures may intuitively suggest a pure phase or a continuous transition for the canonical model. 
This intuition is however false; the characteristics and order of the transition do not depend on the statistical ensemble, 
and the phase diagram of figure 2] is still pertinent to the canonical ensemble 17]. Indeed the relation between the two 
ensembles can be written as 

log PpMt) = log Z () {A t ) + (3pA t - log Zp„ (1) 

where Zp^, Zp(A t ) are the partition sums in the two ensembles. Eq.Q) shows that in the whole region where the 
grancanonical distribution Pp^(At) is convex, the canonical equation of state 

_ 1 d log Z (i 

AW ~ dA t () 
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presents a back bending, which is an unambiguous signal of a first order phase transition 18]. At each temperature 
the maxima of Pp^ correspond to the two ending points of the tangent construction for eq.||2J), i.e. to the borders of 
the coexistence region in the canonical ensemble. 

The qualitative behavior of Am (T) in the canonical ensemble does not change with the density of the system. In 
particular, the Am fluctuation passes systematically through a maximum. The locus of these maxima is displayed on 
the phase diagram in figure ^ We can see that the maximum fluctuation approximately corresponds to the transition 
temperature only at the critical point. At subcritical densities this maximum lies inside the coexistence region of 
the first order phase transition. The results of figure [21 show that the double hump criterium for a first order phase 
transition does not hold if a constraint is put on a variable closely correlated to the order parameter under study. 

IV. CONSERVATION LAWS AND DELTA SCALING 

We can ask the question whether a detailed study of the scaling properties of the Am distribution may give extra 
information on the transition and discriminate first and second order. Following the arguments of ref.|3( we consider 
the distribution 

= $ { % M ^ M ) ^ {AM)AP{AM) (3) 

where A* M is the most probable value of Am and < A < 1 is a real number. At a continuous phase transition point 
the distribution of the order parameter is expected to fulfil the first scaling law, i.e. the distribution $ should be scale 
invariant with A = 1. The scale invariance of <E> for a given value of A is generically refered to as A scaling, and the 
transition observed experimentally[5j from a A w 1/2 to a A« 1 scaling by varying the centrality of the collision and 
therefore the energy deposited in the system, has been taken as a signal of a continuous phase transition. 

A practical difficulty in testing A scaling is that for a given distribution the value of A that corresponds to scale 
invariance, if any, cannot be known a-priori. This difficulty can be circumvented using the fact that the scaling © 
imposed < Am > a oc a. Then it is immediate to verify that eq.(|3J) can be equivalently written as the ensemble of the 
two conditions 

o Am P{A m ) (4) 

K(A M ) 2A (5) 

where ^ is a scale invariant distribution and K is a constant. Since in presence of a scaling, the difference between 
the most probable A* M and the average < Am > scales like a, Am can either be one or the other. In the later case the 
occurence of A scaling study coresponds to the invariance of the centered and reduced distribution. If this distribution 
does not show scale invariance, we can exclude the existence of any A scaling law. If the function \I/ is scale invariant, 
this corresponds to a A scaling if and only if the log- log correlation between the average and the variance is linear; 
in this case the slope of the correlation gives the value of A. The practical advantage of testing eas. (|4l5|l instead of 
eq.@ is that we can check scale invariance without any a-priori knowledge of A. 

The standard way of testing scale invariance is to consider a specific point of the phase diagram and consider the 
centered and reduced Am distributions obtained by varying the size of the lattice and, as a consequence, the total 
number of particles. For the canonical case at the thermodynamical critical point, this analysis is shown on the left 
side of figure |3J Both eq.(QJ and eq.© are well verified, in agreement with the expectation of a first scaling law at a 
continuous transition point |3(. A comparable quality scaling is however observed also at subcritical densities at the 
temperature corresponding to the maximum Am fluctuations (right side of figure |3J) . This finding is in agreement 
with ref.|ll|. Together with the analysis of the phase diagram this means that such a scaling also approximately 
applies in the coexistence region of a first order phase transition, if the order parameter is not free to fluctuate but is 
constrained by a conservation law. 

V. DELTA SCALING AS A FUNCTION OF THE SYSTEM EXCITATION 

In the experimental application to nuclear multifragmentation[5j the system size cannot be varied as freely as in the 
lattice gas, since the maximum size for a nuclear system is of the order of 400 particles. To explore different values of 
(Am), the same system has been studied at different bombarding energies 5] and/or different impact parameters|2l). 
In a similar way, we have kept the total number of particles constant and we have varied the temperature. To fix the 
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FIG. 3: Delta scaling analysis for the canonical lattice gas model at the critical point p = po/2,T — T c (left side) and inside 
coexistence p — po/i at the point of maximal Am fluctuation (right side). Upper part: centered and reduced Am distributions. 
Lower part: correlation between the first two moments and linear interpolation according to eq.J^J. The linear size of the 
lattice is varied as L — 5, 6, 8, 12, 16. 



ideas, we have chosen the simplest thermodynamical path from coexistence to the fluid phase passing through the 
critical point, p(T) = cte = po/2. The resulting <3> functions are displayed in figure 0] No scaling is observed: the 
function 'F continuously evolves from a distribution with a tail extending towards the low mass side compared to the 
average at low temperature while the opposite is true at high temperature. 

If we look at the behavior of the variance as a function of the first moment, the log-log correlation is nowhere linear 
showing that the large fragment fluctuation does not evolve like a power of the average fragment size. The bell shaped 
behavior of this curve is due to the mass conservation constraint, that forces the fluctuation to vanish both at low 
and at high (Am) values. The observed maximum is in fact the maximum fluctuation point shown in figure 1 which, 
at the critical density, occurs close to the critical point and, for sub-critical densities, is located inside the coexistence 
region. 

To qualitatively compare with experimental A-scaling analysis, we have to remember that the studied experimental 
distributions only cover the multifragmentation regime and do not explore the decreasing part of the ua m {{Am)) 
correlation which would correspond in the nuclear case to evaporation from a compound. Focusing now on the 
fragmentation region, we show in figure 0] the best power law interpolations of the average and variance correlation to 
be compared with the published experimental analysis presenting a A = 1 to a A = 1/2 regime. In this interpretation, 
the crossing point between the two power-law fits is interpreted as a "transition" point. By construction it appears 
to be at a higher temperature then the maximum fluctuation which at this critical density comes out to be close to 
the critical temperature. 

To better study the possible occurrence of a power law scaling of the large fragment fluctuation we can study the 
double logarithmic derivative 

A ' = #T (6) 

d{A M ) 

In presence of a A-scaling this quantity should be constant. Figure 21 shows that A' is a smoothly decreasing function 
passing through the values 1 and 1/2 before going through at the maximum fluctuation point and then becoming 
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negative as a consequence of the mass conservation law. No plateaus of A' are observed confirming the absence of 
scaling. 




(A M -<A M >)/a log(<A M > 2 ) 

FIG. 4: Delta scaling analysis for the canonical lattice gas model at constant density p — po/2 varying the system temperature. 
Left part: centered and reduced Am distributions for temperatures varying from T = Lie to T = 1.25e. Right part: correlation 
between the hrst two moments and linear interpolations according to eq.@. Temperatures range from T = .3e to T = 2.3e. 
The vertical lines indicate the temperature of maximum Am fluctuations and the critical temperature. The double logarithmic 
derivative A' is shown in the inserted figure. 

This violation of scaling occurs in spite of the fact that a continuous phase transition point (the thermodynamic 
critical point) is explored in the simulations. At this point the distributions indeed follow the first scaling law (left 
part of figure |3J| but this information is lost if the different distributions are generated by varying the temperature. 
This is not only true for the transition point, but also for the supercritical regime. Indeed this regime has been shown 
to exhibit the second scaling law A = 1/2 in the Potts modelQ (or something close to it A « 0.6 for the IMFM[Ti|L 
while in the representation of figure 01 scaling can everywhere be excluded. 

The conclusion is that scale invariance can only be tested by varying the total system size. However other information 
on the phase transition can be accessed through the study of the Am distribution with a fixed total number of particles, 
as we now show. 



VI. SIGNALS OF PHASE TRANSITION AND OF ITS ORDER 

The first two moments of the distribution in the canonical ensemble and the corresponding most probable value 
A* M are displayed in figure [S] for three different densities. Let us look at the p < p c case first. If the first and second 
moment show smooth behaviors dominated by the conservation law constraint, the transition is still apparent in the 
behavior of A* M which rapidly changes at a temperature close to the transition point. This sudden decrease is due to 
a change of sign in the asymmetry of the distribution. As such, the qualitative behavior of A* M (T) is independent of 
the density. A great number of continuous transition signals has been observed in different mass conserving models at 
densities that do not correspond to a continuous phase transition |ri PI Isl Wi TToL lLl| . The same happens for the most 
probable value of Am- This variable shows for all densities a sudden drop at a temperature T t (p) which corresponds 
to the maximum of the Am fluctuations. As we have already stressed, these temperatures approximately coincide 
with the transition temperature only at the critical density (see figure ^| . The behavior at supercritical densities 
reflects a geometric phase transition which has no thermodynamic counterpart, while if fragmentation takes place at 
low density the Am drop can be taken as a signal of phase coexistence. 

In order to discriminate between the different density regimes and recognize the order of the phase transition, we 
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FIG. 5: First moments of the Am /At distribution as a function of temperature for a 8x8x8 lattice in the canonical ensemble at 
p/po — I/A (left), p/po = 1/2 (medium), and p/po = 3/4 (right). Upper part: mean value (full line) and variance (dashed line). 
Lower part: most probable value of the Am distribution normalized to the mean. The vertical lines indicate the temperature 
of maximal Am fluctuations and the transition temperature for each density. 



have to quantify the Am fluctuation peak. In the grancanonical ensemble, the A t fluctuation is directly linked to the 
susceptibility via 

X=^ = « t ) 2 (7) 

To work out a similar expression for the canonical ensemble, let us assume that Am and the other fragments are 
statistically independent, i.e. the total density of states is factorized 

W t {A m , A m , E M , E m ) = W M (A M , E M ) ■ W m {A m , E m ) (8) 

where we have defined the total number of particles not belonging to the largest fragment as A m = A t — Am, and 
the corresponding energy E m — E t — Em- This hypothesis is reasonably well verified in the Lattice Gas model, since 
the correlation coefficient between Am and A m in the grancanonical ensemble comes out to be close to zero except in 
the very dense regime p/po ~ 1. The factorization of the state densities implies a convolution of the corresponding 
canonical partition sums 

E A 

Zp(A t ) = [ dE t e~ f3Et [ 1 dE m [ ' dA m ■ W m (E m , A m )W M (E t - E m ,A t - A m ) 
J Jo Jo 

dA m Zf{A M )Zf{A t -A M ) (9) 

where ZL,i = m,M describe the contribution of the largest fragment and of all the others respectively. 
The distribution of the largest fragment reads 

P^aMm) - Zp\A t )Zf{A M )ZJ{A t - A M ) (10) 

A Gaussian approximation of this distribution leads to|l9j| 

^ = ( Xm (A t -A* M ) ~ HM^) 
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where a\ is the fluctuation of the Am distribution and the partial susceptibilities are defined as Xi 1 = ^^iAt)- 
The above derivation is valid for a system whose state density depends on the two extensive variables, number of 

particles A and energy E. In the case of the fragmentation transition a third extensive variable, the volume V, has 

also to be considered. We show in the appendix that in this more general case eq. Ijllf) can still be derived, but a 

dilute limit Vt = V m + Vm ~ V m has to be considered. 

According to the general definition of phase transitions in finite systems 0, Hjj, the generalized susceptibility 

associated to an order parameter is negative in a first order phase transition in the statistical ensemble where the 

order parameter is subject to a conservation law. We therefore expect a negative \m at subcritical densities. Imposing 

Xm < in ea. ((TTJl leads to 



Comparing to eq. J7J) this finally gives 



> /3 l Xm{A t -A* M ) 



) 2 



(12) 



(13) 



Equation (|13|) associates the first order phase transition in the canonical ensemble to "abnormal" Am fluctuations, 
in the same way as abnormal partial energy fluctuations sign a first order phase transition in the microcanonical 
ensemble|l9j. The canonical and grancanonical fluctuations are compared in figure for three different densities. 
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FIG. 6: Ratio between the grancanonical and canonical fluctuation of the number of particles A m not belonging to the largest 
cluster, as a function of the temperature for a 8x8x8 lattice at p/po = 1/4 (left),p/po = 1/2 (medium), and p/ po = 3/4 (right). 
Dashed lines: average canonical Am values normalized to the total number of particles. Vertical lines in the two left panels: 
limit of the region of negative susceptibility from the canonical p,(A) equation of state eq. |^J) ||. 



Independent of the system density our approximation ea. l|13|) turns out to be incorrect at very low temperatures, 
when the average size of the largest cluster (dashed lines) exceeds about 80% of the available mass. In this case the 
hypothesis of statistical independence between A m and Am cannot be justified and the canonical mass conservation 
costraint trivially reduces the canonical fluctuation. However as soon as the average Am value drops, we can see that 
the region of negative susceptibility can be well reconstructed through ea.l |13(l . and in particular its border (vertical 
lines) is very precisely determined by the equality condition between the two fluctuations. At supercritical densities the 
dilute gas approximation we have employed breaks down independent of the temperature and the susceptibility cannot 
quantitatively be estimated from the fluctuation signal, however in this regime the relative fluctuation observable does 
not present any peak while only inside the spinodal region of the first order phase transition the canonical fluctuation 
exceeds the grancanonical one. It is clear that this observable allows a unambiguous discrimination between the 
supercritical regime and phase coexistence. 
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VII. CONCLUSIONS 



To conclude, in this paper we have discussed the role of the largest fragment in the framework of the lattice gas 
model. We have shown that this variable can be taken as an order parameter of the frag mentation phase transition if 
this latter belongs to the liquid gas universality class. It has been already observed [ll|,|22] that the phase transition 
can be tracked from the sudden drop of Am close to the transition temperature. This drop is well fitted by a power 
law with a exponent close to the expected value for the liquid gas universality class but finite size effects blur 
the behavior considerably for system sizes comparable to accessible nuclear sizes. However, when no constraints are 
affecting the fluctuations of the order parameter such as in the grand canonical ensemble, we have shown that the 
transition is very well defined if instead of the average we look at the most probable value of Am- Indeed crossing a 
first order phase transition point this variable is discontinuous independent of the system size. The important result 
is that if we look at this variable finite size effects do not constitute a major problem to identify a phase transition 
nor to recognize its order. 

On the other hand important ambiguities arise from the non equivalence of statistical ensembles inside a phase 
transition. Indeed the distribution of the order parameter is strongly deformed by the presence of conservation laws 
in the system under study. If we look at Am as an order parameter, the double hump criterium for a first order phase 
transition does not apply any more in the canonical or microcanonical ensemble because of the strong correlation 
between the conserved total number of particles and the order parameter. The mass conservation constraint induces 
a maximum in the fluctuation of Am which is not necessarily correlated with the properties of the phase diagram. 
We observe maxima both at the critical density close to the critical point and at sub-critical densities inside the 
coexistence zone. Moreover, the presence of this maximum can simulate a transition from a A = 1 to a A = 1/2 
scaling law in a region above the maximum fluctuation. It is clear that other observables have to be employed if we 
want to conclude about the order and nature of the phase transition. One such observable is the numerical value of 
the fluctuation of Am, which is by construction identical to the fluctuation of the number of particles that do not 
belong to the largest cluster A m : if and only if the system crosses the phase coexistence region of a first order phase 
transition, this fluctuation overcomes the corresponding value in the grancanonical ensemble. 



VIII. APPENDIX: DERIVATION OF EQ.(IO) 

The density of states is a function of all the relevant extensive variables of the system. For the lattice gas model 
this means W = W(E, A, V). If the largest fragment Am is statistically independent from the other clusters then 

W t (A M , E M , V M ,A m , E m , V m ) = W M (A M , E M , V M ) • W m {A m , E m , V m ) (14) 

where we have defined the total number of particles not belonging to the largest fragment as A m = A t — Am, and 
the corresponding energy and volume E m = E t — EM,V rn = Vt — Vm- Let us first consider the case of an external 
temperature T = /3 _1 and pressure p = j3X. Using the standard definition of the canonical isobar partition sum 

Z/3x(A) = J dEe- pE J dVe- pV W(E,A, V) 

the total partition sum can be written as 

Zpx{A t ) = [ dEte-*** [dV t e- pVt dE m dV m [ ' dA m • 
J J Jo Jo Jo 

W m (E m ,A m , V m )W M (E t -En, At- A m , V t - V m ) 

or equivalently 

Z (A t ) = dA m Zf{A M )Z^{A t - A M ) (15) 



where Zp, i = m,M describe the contribution of the largest fragment and of all the others respectively. In the isochore 
case V m + V)vf=cte, the convolution of the partition sum is less straightforward because of the presence of the volume 
integral 



A t r V t 

M 



Z (A u V t ) = / dA m / dV m Z^ (A M , V M )Zf{A t - A M , V - V M ) (16) 
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Let us introduce the partial pressures p, = (5~ x 9l gy f (A*, V*) and chemical potentials y>i = j3~ x ^q^ (A'i K*) a ^ the 
most probable volume and mass partition A* M , V M . Equilibrium between the two components implies [i m = /im , 
Pm — Pm- A saddle point approximation then gives 

Zf ZJ » exp (-/3 L4^/ M + A* m f m }) ■ 



exp 
exp 



\ (A m - A* m f (XM + Xm ) + \ (Vm - V„*f (Km 1 + K- 1 ) ) 

( dp M , dp n 



A* m ) (V m ~ V*) 



\dA 



M 



dA n 



where fa — —TlnZl(A*)/A*,i — m,M are the most probable free energies per particle, the partial susceptibilities 

and compressibilities are defined as \i X — §a~(A* ,V*), k^ 1 — §y-(A* ,V*), and the conservation constraints make 
the linear terms vanish. In the dilute limit V t = V m + Vm ~ V m the density variation of the "gas" component m is due 
to its number variation dp m = dA m /Vt and the volume variation can be neglected respect to the number variation 
V m - V m < A m - A* m giving 



Z p (A u V t ) « / dA m Zf(A M , C)4 m) {At - A M , V t - V*) 
Jo 

Both in the isobar (|15|) and in the isochore (|17f) case, the distribution of the largest fragment reads 

P P A t {A M ) = Z- 1 Z^(A M )Z^(A t - A M ) 
Implementing the saddle point approximation we can identify 

1 1 



Xm(A t -A* M ) xm{A* m ) 



(17) 



(18) 



(19) 



where <j\ is the fluctuation of the Am distribution. 



X. Campi, J.Desbois, E.Lipparini, Phys. Lett. 138B (1984) 353. 
K. Binder, D. P. Landau, Phys. Rev. B 30 (1984) 1477. 
R. Botet, M. Ploszajczak, Phys. Rev. E62 (2000) 1825. 

B. Tamain et al., Nucl.Phys.A, in press; M.F.Rivet et al, nucl-ex/0412007 
R.Botet et al., Phys.Rev.Lett. 86(2001)3514. 

C. B. Das, S. Das Gupta and A. Majumder, Phys. Rev. C65 (2002) 34608; C. B. Das et al., Phys. Rev. C66 (2002) 044602. 
K. A. Bugaev et al., Phys. Rev. C62 (2000) 044320 and Phys. Lett. B 498 (2001) 144. 

F.Gulminelli, Ph.Chomaz, Phys.Rev.Lett.82 (1999) 1402; Ph.Chomaz, F.Gulminelli, Int. Journ. Mod. Phys. E8 (1999) 527. 

F. Gulminelli et al., Phys. Rev.C (2002) 51601. 

A. H. Raduta et al., Phys. Rev. C 65, 034606 (2002). 

J. M. Carmona, J. Richert, P. Wagner, Phys. Lett. B531 (2002) 71. 

C. N.Yang, Phys.Rev.85 (1952)809 

A. Coniglio and W. Klein, J. Phys. A13 (1980) 2775; X. Campi, H. Krivine and A. Puente, Physica A 262 (1999) 328. 
F. Gulminelli and Ph. Chomaz, Phys.Rev.E 66 (2002) 046108. 

K. C. Lee, Phys. Rev. 53 E (1996) 6558; Ph. Chomaz, F. Gulminelli, Physica A (2003). 
Ph.Chomaz, F.Gulminelli, V.Duflot, Phys.Rev.E64 (2001) 046114. 
F. Gulminelli et al., Phys. Rev. E 68 (2003) 026120. 

D. H. E. Gross, " Microcanonical thermodynamics: phase transitions in finite systems", Lecture notes in Physics vol. 66, 
World Scientific (2001). 

P.Chomaz and F.Gulminelli, Nucl. Phys. A647 (1999) 153. 

Ph.Chomaz and F.Gulminelli, in 'Dynamics and Thermodynamics of systems with long range interactions', Lecture Notes 

in Physics vol.602, Springer (2002). 

J.Frankland et al., nucl-ex/0404024 unpublished. 

M. Pleimling and A. Hueiler, J. Stat. Phys. 104 (2001) 971. 



